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Abstract 

We use a recently improved version of the chiral nucleon-nucleon potential at next-to-next-to- 
leading order to calculate the pairing gap in isospin-symmetric nuclear matter. The pairing 
potential consists of the long-range one- and two-pion exchange terms and two short-distance 
NN-contact couplings. We find that the inclusion of the two-pion exchange at next-to-next-to- 
leading order reduces substantially the cut-off dependence of the ^Sq pairing gap determined by 
solving a regularised BCS equation. Our results are close to those obtained with the universal 
low- momentum nucleon-nucleon potential Viow-k or the phenomenological Gogny DIS force. 

PACS numbers: 13.75.Gx, 21.30.Cb, 21.30.Fe, 21.60.-n 
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The self-consistent mean-field framework, extended to take into account the most impor- 
tant correlations, provides at present the only viable microscopic description of structure 
phenomena in light and heavy nuclei over the entire periodic table. A broad range of suc- 
cessful applications to nuclear structure and low-energy dynamics characterizes mean-field 
models based on the Gogny interaction^the Skyrme energy functional, and the relativistic 
meson-exchange effective Lagrangian The effective forces used in these models contain 

a moderate number of free parameters that are adjusted to global properties of a small set 
of spherical and stable nuclei, rather than to the observables of free NN-scattering. In other 
words, even though the global effective nuclear interactions model the interaction between 
nucleons in the nuclear medium, they are not necessarily related to any given NN potential. 

On the other hand, a completely new framework has recently been developed in which nu- 
clear interactions are formulated in terms of effective field theories 3,0,0]. The key element 
is the separation of long- and short-distance dynamics and an ordering scheme in powers of 
small momenta. The nucleon-nucleon potential, constructed in the chiral perturbation the- 
ory, consists of the long-range contributions generated by one-, two- and three-pion exchange 
[a, I2I , and a set of contact-terms encoding the short- distance NN-dynamics. The associated 
low-energy constants are adjusted to the empirical NN-phase shifts and deuteron properties. 
It has been shown that their values can be understood in terms of the heavy-mass reso- 
nance exchanges jsj. Furthermore, when the isospin-breaking corrections are systematically 
included, the chiral nucleon-nucleon potential reaches almost the same accuracy as the more 
phenomenological "high-precision" nucleon-nucleon potentials. This has been demonstrated 
in numerous calculations of two- and few nucleon systems. 

Another line of approach, developed by the Stony-Brook group, applies renormalisation 
group arguments to eliminate the high-momentum components from phase-shift equivalent 
nucleon-nucleon potentials. By integrating out the high-momentum components below a 
scale A ^ 2 fm~^, a "universal" low- momentum nucleon-nucleon Viow-k potential emerges j^. 
Consequently, this potential operates only in the subspace of nucleon states with momenta 
p < A. The pairing properties of nuclear matter derived from the Viow-k potential have 
been investigated recently in Ref. Q]. Good agreement with the phenomenological Gogny 
pairing interaction has been found for a wide range of nuclear densities. Calculations of 
nuclear matter properties have also been performed using the Viow-k potential. The results 
obtained so far in the Hartree-Fock and Brueckner-Hartree-Fock approximations are however 
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unsatisfactory because no saturation occurs in the nuclear matter equation of state 111]. It 
has been concluded that the Brueckner-Hartree-Fock approximation, based on the l^ow-k 
potential, is applicable only at very low densities. These findings are consistent with the 
fact that nuclear pairing is known to be primarily a low-density phenomenon. The maximum 
of the ^5*0 pairing gap typically lies at a density p ^ pq/A, where po denotes the nuclear matter 
saturation density. The purpose of the present paper is to investigate the pairing properties 
of the chiral nucleon-nucleon potential. 

Quite generally, the momentum and density-dependent pairing field A{k,kf) in infinite 
nuclear matter is determined by the solution of the BCS gap equation 

*" ^° ^J\E(p,k,)-£(k,,k,W + A(p,k,r- 

V{p, k) represents the off-shell pairing potential in momentum space, £{p, kf) is the quasi- 
particle energy, and £{kf, kf) is the Fermi energy. 

The effective force in the pairing channel is generated by the sum of all particle-particle ir- 
reducible Feynman diagrams 3]- In niost application to nuclear and neutron matter, 
however, only the lowest-order term, which corresponds to the bare nucleon-nucleon inter- 
action, is retained Terms of higher order in the effective pairing interaction represent 
screening corrections to the bare force, caused by medium polarization effects. Numerous 
studies have shown that polarization effects can have a pronounced influence on the cal- 
culated values of the pairing gaps (see, for instance, Refs. |l8|). The influence of 
both the vertex corrections to the pairing interaction, and the self-energy corrections, on 
the properties of ^5*0 pairing in neutron and nuclear matter has recently been studied in the 
framework of the generalized gap equation 2^. It has been found that the two effects 
lead to a strong suppression of the pairing correlations in neutron matter (the pairing gap 
is reduced by more than 50% with respect to the BCS prediction), whereas they cancel each 
other out to a large extent in symmetric nuclear matter. The pairing correlation energy in 
finite nuclei can be calculated in the local density approximation (LDA). LDA calculations 
with the Gogny Dl force have been compared with exact Hartree-Fock-Bogoliubov (HFB) 
calculations of the pairing correlation energy for many spherical nuclei Except for shell 
effects, the results of these calculations are in close agreement, and thus one should not 
expect large effects from medium corrections not included in the BCS limit. In the present 
analysis we only consider symmetric nuclear matter and, therefore, medium polarization 



effects will not be taken into account at this stage. 

For the pairing potential V{p, k) we employ the improved version of the chiral nucleon- 



nucleon potential derived in Refs. j22l . |23| | . An overly strong medium-range attraction in the 
isoscalar central part of the chiral two-pion exchange at the next-to-next-to-leading order in 
the chiral expansion, present in earlier versions of the chiral potential P], has been removed 
by using the spectral function regularization method. Essentially, this means that only 
TTTT-intermediate states of invariant mass below a scale A (where the chiral effective field 
theory is applicable) are taken into account in the pion-loop integrals, while shorter range 
contributions should be represented by NN-contact interactions. 

The bare chiral NN-potential used in the present study consists of the one-pion exchange 
term, the NN-contact interaction, and the irreducible two-pion exchange term 

V{p, k) = V^^""^ {p, k) + V'^"'^ {p, k) + 1^(2") (p, k) . (2) 

This approximation for the two-nucleon potential is valid for small values of the momenta 
p and k and it breaks down for momenta above the chiral symmetry breaking scale. An 
additional cut-off A, included by multiplying the potential V{p, k) with a regulator function 

V{p,k)-.f\p)V{p,k)f\k), (3) 

prevents the growth of the potential with increasing momenta p and k. Following the 
procedure of Ref. we employ the exponential regulator function 

f\p) = exp(-//A') • (4) 

Even though both cut-off parameters: A and A, are introduced to remove high momentum 
components of the interacting nucleon and pion fields, their roles are different. The inclusion 
of A removes the short-distance portion of the TPE component of the nuclear force, whereas 
the parameter A ensures that high-momentum nucleon states do not contribute to the scat- 
tering process. If all terms in the EFT expansion are included, low-energy observables 
should not depend on the cut-off parameters. In practice, however, the expansion is always 
truncated at some order. Consequently, the observables depend on the cut-off parameters 
to some extent, but this dependence becomes weaker when higher-order terms are included. 
Unless stated otherwise, we use the value A = 550 MeV in the following calculations. 
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The one-pion exchange contribution in the S-wave channel reads 

The following numerical values have been used for the nucleon axial vector coupling constant 
qa, the weak pion decay constant and the pion mass m^: qa = 1-3, = 92.4 MeV, 
= 135 MeV. 

The NN-contact interactions includes the contact operators without and with two deriva- 
tives 

K(-)(p,A;) = ^ + ^(/ + fc^). (6) 
Ztx Ztt 

For each choice of the cut-off parameters A and A, the low-energy constants Cisg and Cisg 
are determined by a fit to the ^5*0 phase shifts below the inelastic A^A^vr-threshold. The cut- 
off dependence entering through the regulator function /a(p) is to a large extent cancelled 
by that of the running low-energy constants Cis^, and Cis^- In the present calculation we 
employ the central values of the low-energy constants at the next-to-leading order (NLO) 
and next-to-next-to- leading order (NNLO) listed in Table 4 of Ref. Q]. These particular 
values of the low-energy constants have been adjusted for A = 600 MeV and A = 550 MeV. 

The irreducible two-pion exchange contribution is most conveniently calculated via a 
twice subtracted dispersion relation 

y(2-)(p,A;) = - dfilm[Vc + Wc-2fi^VT-2ij,^WT] 

I /i /i^ pk fi'^ + {p — kY J 

ImVc, ImyTc; ImVr, and ImPVr are the imaginary parts of the isoscalar central, isovector 
central, isoscalar tensor, and isovector tensor NN-amplitudes. The weighting function in 
curly brackets originates from projecting the NN-potential onto the spin-singlet S-wave. 
The cut-off parameter A, which restricts the spectral functions calculated in the chiral 
perturbation theory to their low-energy domain of validity, is kept fixed throughout this 
work at the value A = 600 MeV. The imaginary parts of the isovector central and isoscalar 
tensor NN-amplitudes contribute at the NLO, while the imaginary parts of the isoscalar 
central and isovector tensor NN-amplitudes contribute at the NNLO in the chiral expansion. 
Explicitly, the spectral functions corresponding to the two-pion exchange at the NLO read 



37r/x(4A)^ 



4m^(l + - hg\) + ii\2?>g\ - IQg^ - 1) + 
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and at the NNLO 
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ImVcitfi) = (4ci - 2c3)m^ + C3 /i' (2m^ - /i^) , (10) 



^2 



WTM = ^^(4m^-/i^), (11) 
with the low-energy constants: Ci = — 0.81GeV~^, C3 = — 3.40GeV~^, and C4 = 3.40 GeV~^ 
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For the single-particle spectrum which enters the gap equation (p^, we employ the simple 
quadratic form 

S(p,kf)-£(kf,kf) = . (12) 

v/^' }) V /; 2M*{kf) ^ ' 

Such an approximation is sufficient since momenta p around kf give the dominate contri- 
bution to the integral in Eq. (^. The effective nucleon mass M*{kf) is deduced from a 
very recent three-loop calculation of nuclear matter in chiral perturbation theory including 
explicit A(1232)-isobar degrees of freedom Q]. 

As we have already emphasized, the low-energy observables should not depend sensitively 
on the precise value of the cut-off parameter A. In addition, the dependence should become 
weaker when higher order terms are included in the small momentum expansion. To verify 
this we have calculated the pairing gaps Aj = A{kf,kf) in symmetric nuclear matter for 
several values of the cut-off parameter A between 450 MeV and 650 MeV. For the low-energy 
constants Cis^ and Cis^ we have used the central values listed in Table 4 of Ref. [2^, and 
the spectral function cut-off A has been fixed at A = 600 MeV. Notice that no attempt 
has been made to readjust the LECs when changing the cut-off parameter A. The results 
are displayed in Fig. ^ At the NLO the pairing gap strongly depends on the value of A. 
With increasing the cut-off from A = 450 MeV to A = 600 MeV, the maximum value of 
the pairing gap increases from A™"^ = 0.7 MeV to Aj^""^ = 5.81 MeV, and the value of the 
Fermi momentum for which the pairing gap reaches maximum changes from pj"-^ = 0.65 
fm~^ to pY"-^ = 0.95 fm~^. For A = 650 MeV the maximum value of the gap is well above 10 
MeV and the calculation obviously diverges. The results are much improved at the NNLO: 
the position of the maximum changes very little with increasing A, and the value of the 
pairing gap at maximum increases less than 0.8 MeV between A = 450 MeV and A = 650 
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MeV. Since the gap at the Fermi surface is very sensitive to changes in the effective pairing 
potential, this indicates rapid convergence for the EFT expansion. 

In the left panel of Fig. |21 we compare the density dependence of the pairing gaps in 
nuclear matter, calculated by employing the chiral NN-potential and with the effective Gogny 
interaction. The solid and the dashed curves correspond to the chiral potential at the 



NLO and the NNLO, respective 
with the DIS Gogny interaction 



The dash-dotted curve is the pairing gap calculated 
26| . Similar to the analysis of Ref. Ql? where the Gogny 
pairing gaps were compared with those calculated with the low-momentum V^ow-fc interaction 
extracted from realistic forces, we have calculated the pairing gaps using both the non- 
interacting single-particle spectra (heavy curves), and the medium- modified single-particle 
spectra (light curves). The medium- modified single-particle spectra are computed in the 
Hartree-Fock approximation for the Gogny interaction, and by using Eq. (fT^ for the chiral 
NN-potential. First we notice that for the noninteracting single-particle spectra there is a 
qualitative and, at low densities kf < 0.8 fm~^, even a quantitative agreement between the 
pairing gaps calculated with the chiral NN potential and the effective Gogny interaction. 
The inclusion of the NNLO contribution increases the value of the pairing gap at maximum 
and shifts the position of the maximum to a slightly higher value of kf, in closer agreement 
with the Gogny pairing gap. For kf > 1 fm~^ and closer to saturation density, we find a 
more pronounced difference between the Gogny pairing gap and the gaps calculated with 
the chiral NN potential. This difference with respect to the Gogn y p airing gap at higher 
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densities has also been observed for other bare NN interactions 

The effective nucleon mass is reduced in the nuclear medium. This results in a lower 
density of states around the Fermi surface and, hence, in weaker pairing correlations. The 
reduction of the pairing gap is much more pronounced for the Gogny effective interaction 
than for the chiral NN potential. This is caused by a much stronger reduction of the effective 
mass for the former interaction in the relevant range of densities 0.4 < kf < 1.2 fm~^ The 

" n 

effective nucleon mass M*{kf) / M from the chiral EFT three-loop calculation [2^, and for 
the DIS Gogny effective interaction are shown in the right panel of Fig. |21 The density 
dependence is very different in the two cases, and this means that one should not, at least 
on the quantitative level, compare the corresponding pairing gaps calculated with medium- 
modified single-nucleon spectra. 

In Fig. El we plot the momentum dependence of the pairing potentials V{p, k) for three 
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different values of the momenta k = 0.4,0.8,1.2 fm~^. While the Gogny interaction and 
the chiral NN-potential at the NLO display only a qualitatively similar momentum depen- 
dence, the inclusion of the NNLO contributions brings the chiral NN-potential in remarkable 
agreement with the Gogny DIS interaction. The same is true and even more pronounced, 
of course, for the diagonal matrix elements V{k, k) shown in the right panel of Fig. |21 A 
similar result has been obtained for the V^ow-k potential in Ref. jlo| . 

In conclusion, we have presented a study of the chiral nucleon-nucleon potential pairing 
properties. The potential consists of the long-range one- and two-pion exchange terms, and 
it includes the short-distance NN-contact terms. Our main result is that the inclusion of 
the two-pion exchange at the next-to-next-to-leading order reduces substantially the cut-off 
dependence of the ^5*0 pairing gap, indicating good convergence of the small momentum ex- 
pansion. In addition, the inclusion of the NNLO contributions brings the chiral NN-potential 
in remarkable agreement with the Gogny DIS interaction and the V^ow-k potential. In fu- 
ture we plan to include the NNNLO terms in the chiral expansion of the bare NN-potential 
and extend our study to the neutron matter pairing properties. This will necessitate the 
inclusion of medium modification effects, such as the polarization effects, vertex corrections 
and Pauli-blocking effects. The chiral effective field theory at finite density represents a 
systematic framework for such extension. Work along these lines is in progress. 
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FIG. 1: Pairing gap at the Fermi surface Af = A(A;/, kf) as function of the Fermi momentum, for 
different values of the cut-off parameter A, at NLO (upper panel) and NNLO (lower panel). 
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FIG. 2: Pairing gaps at the Fermi surface = A{kf, kf) in symmetric nuclear matter, calculated 
with the chiral NN-potential @ at NLO and NNLO, and with the effective Gogny interaction 
DIS 26]. The heavy and light curves refer to calculations with non-interacting and medium- 
modified single-particle spectra, respectively. The panel on the right illustrates the density de- 
)endence of the effective nucleon masses for the Gogny interaction and for the chiral potential 
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FIG. 3: Momentum dependence of the pairing potentials V{p, k): chiral NN-potential at NLO and 
NNLO, and the Gogny DSl interaction, for three different values of the momenta A; = 0.4,0.8, 1.2 

fm^^. In the panel on the right the corresponding diagonal matrix elements V{k,k) are shown as 
functions of the momentum k. 
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